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Abstract 

We present a geometric, model- independent, argument that explains why the Tsallis en- 
K^ tropy describes dynamical systems having vanishing largest Lyapunov exponent. We 

t^^ ' employ the Jacobi/geodesic deviation equation for an effective negative curvature Rie- 

O ; 

t^^ I mannian metric reflecting the Tsallis entropy composition property, whose solution gives 



the desired result. Extending the essential parts of the argument from Riemannian man- 
O ■ ifolds to CAT(k), k < spaces, we see that the conclusion remains valid in the case 

of interacting systems described by different entropic parameters. This conclusion is in 
agreement with all currently known results. 
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1. Introduction 

The Tsallis entropy is a single-parameter family of functionals, first introduced in the 
Physics literature in 1988 [1], providing an alternative to the Boltzmann/Gibbs/Shannon 
(BGS) entropy used in the statistical description of a system. Consider a probability 
distribution {pi}, i G / in a discrete sample space indexed by / C N. Its Tsallis entropy 
is defined as 

It may be worth comparing (1) with the BGS entropy 

SBGsiiPi}) = -ke^Pi log Pi (2) 

where we immediately observe as the non-extensive/entropic parameter g — )> 1, we get 

lim Sg = Sbgs (3) 

An analogous definition of the Tsallis entropy can be given for continuous sample spaces. 
Henceforth, we will be setting the Boltzmann constant ks = 1, for simplicity. 

The Tsallis entropy, conjecturally, describes collective phenomena with long-range spa- 
tial and temporal correlations [2], [3], systems whose phase portraits exhibit a fractal-like 
behavior etc. for which there is no reason or justification why their description by the 
BGS entropy should be accurate or even valid [3] (and references therein). Following the 
approach and viewpoint of Boltzmann [4], [5], one can state that the dynamical basis 
of the Tsallis entropy remains unclear so far [3]. This mirrors the existing difficulties in 
deriving the BGS entropy [3] - [6] from dynamical principles. As in the latter case, in the 
case of the Tsallis entropy, some progress has been made in identifying characteristics of 
systems effectively described by it. One class of such systems are ones exhibiting "weak 
chaos" or being at the "edge of chaos" for some subset of their parameter space [2], [3], 
[7] - [10]. This is more accurately expressed by referring to them as dynamical systems 
having vanishing largest Lyapunov exponent [3], [7] - [10]. 

A general, model- independent and formal justification of why such systems are de- 
scribed by the Tsallis entropy has been lacking so far [3], [11]. This is the issue that 
we are addressing in the present work. To be more specific, we present an argument 
about the converse: we argue, in Section 2, that if a system is described by the Tsallis 



entropy, then its largest Lyapunov exponent vanishes. This statement is in agreement 
with all currently known resultcl about the underlying dynamics of systems described by 
the Tsallis entropy. In Section 3, we extend this conclusion to interacting systems having 
different values of q. Such interacting systems should be modelled on CAT(k), k < 
spaces, as we pointed out in [12]. Section 4, contains some general comments and points 
toward further implications of the Tsallis entropy composition property (6), employing 
the underlying concept of hyperbolicity. 

For brevity, we do not provide the required background in Riemannian geometry or 
the metric geometry of CAT(k), k < spaces, referring instead to some of the excellent 
references such as [13], [14] and [15], [16] respectively, on these topics. 



2. Riemannian spaces: vanishing largest Lyapunov exponent 

One of the features distinguishing the Tsallis (1) from the BGS (2) entropy is their 
different composition properties. This distinction has profound consequences for the def- 
initions of independence and additivity/extensivity that pervade Statistical Mechanics 
and Thermodynamics [3] (and references therein). Two systems A and B are conven- 
tionally defined to be "independent" if in their statistical description, the corresponding 
probability distributions pa_, Pb have as composition law the ordinary multiplication 

Pa+b =Pa-Pb (4) 

Here and henceforth A + B indicates the compound system formed by combining A 
and B. For such independent systems (4), the BGS entropy is additive 

Sbgs{A + B) = Sbgs{A) + Sbgs{B) (5) 

as can be immediately seen from (2). By contrast, under the same definition of indepen- 
dence (4), the Tsallis entropy (1) obeys the generalized additivity/composition property 

S,iA + B) = S,{A) + SgiB) + (1 - q)S,{A)S,{B) (6) 

This composition motivated the introduction of a generalized addition [17], [18] by 

x®gy = X + y + {1 - q)xy (7) 



*We are grateful to Professor Constantino Tsallis for pointing this out to us. 



What was lacking for sometime though, was an appropriate generahzed counterpart to 
the usual multiplication that would be distributive with respect to (7). Despite some 
initial misses, two such generalized multiplications were introduced independently in [19], 
[20]. Although conjecturally independent, an explicit equivalence between them is still 
lacking. In these works a field isomorphism was, essentially, introduced that was indicated 
by Tq in [20] , where some if its rudimentary metric and measure-theoretical properties 
on M were examined. In our subsequent work [12], we developed further some of the 
metric consequences of r^, < g < 1, initially for M". Motivated by an analogy with 
the translation invariance of the Euclidean metric, we constructed in [12] a Riemannian 
metric g, induced by (6), with components 



and corresponding line element 



where 



1 

e^^te 



ds' = dx' + e~''''dy' (9) 



t = log(2 - q) (10) 



The BGS composition property ("ordinary addition") is expressed, in this formalism, via 
the effective Euclidean metric 

whose corresponding line element is 

ds% = dx^ + dy"^ (12) 

The two metrics (8), (11) were used in [12] to encode and compare, in metric terms, the 
composition properties of the BGS and the Tsallis entropies. 

The metric tensor (8) turned out [12] to have a constant negative sectional curvature 

k = -[hg{2 ~ q)]' (13) 

which provided a geometric interpretation of the non-extensive parameter q G [0,1). 
Endowing the plane M^ with (8), turned it into a re-scaled version of the hyperbolic 
plane H^. Then, we concluded in [12], that due to (6) which gave rise to (8), the Tsallis 



entropy can be thought as a "hyperbohc analogue" of the BGS entropy. 

Without adding any complexity to the subsequent argument, we will consider in the 
formalism instead of just M^, a general Riemannian manifold (M, g) with tangent bundle 
indicated by TM. Let the Levi-Civita connection compatible with g be indicated by 
V. Such a connection is expressed in terms of g by the Koszul formula 

2g(VxF, Z) = -g{X, [Y, Z]) + g(Z, [X, Y]) + g(y, [Z, X]) 

+X[g(F, Z)] - Z[giX, Y)] + y [g(X, Z)] (14) 

where X,Y,Ze TM. The geodesic equation is [13], [14] 

VxX = (15) 

for X G TM tangent to the geodesic. If J G TM indicates a Jacobi field then it 
satisfies the Jacobi/geodesic deviation equation [13], [14] 

Vx^xJ + R{J,X)X = (16) 

Here R{X,Y)Z indicates the Riemann tensor which is defined by [13], [14] 

R{X,Y)Z = VxVyZ -VyVxZ -V[x,Y]Z (17) 

It is extremely difficult to explicitly solve either the geodesic (15) or the Jacobi equations 
(16), except in a few particularly simple cases, our case of interest being one of them. Let 
61,62 be orthonormal vectors, with respect to g, spanning a 2-dimensional subspace 
of TM in a neighborhood of x G M. The sectional curvature of M in this 2-plane 
subspace of TM is defined by [13], [14] 

A; = g(i?(6i,62)62,6i) (18) 

In our case of interest, initially at least. A; < is constant. Then the Jacobi equation 
(16) reduces to the ordinary differential equation 



which has the general solution 

n— 1 



kJ{s) = (19) 



J{s) = y^ \ ttj exp{\/^-k s) + bj exp(— -\/^ s) \ 6j(s) (20) 



where a^, bi are constants, {ei{s)} are parallel orthonormal vectors (Fermi basis) and 
s is the arc-length parameter of the geodesic whose tangent is X{s). The summation 
takes place over the n — 1 directions orthogonal to X{s). Substituting (10) into (20), 
we find that 



n-l 



J{s) = 2^{a^{2 - q)e' - 6.(2 - g)e^} e,(s) (21) 

Hence, we see that in a Riemannian manifold of constant negative curvature (13), the 
nearby geodesies deviate from each other exponentially in terms of the geodesic arc-length 
s or any of its affine re-parametrizations. The Riemannian metric (8) is a special case of 
the general g, which is defined on M = M^. This exponential deviation of the nearby 
geodesies of a manifold of negative sectional curvature should be contrasted to those of 
the Euclidean metric (11) which has k = 0: in the Euclidean case the geodesies separate 
linearly as functions of the arc-length parameter s, as can be immediately seen from 
(19) by setting k = 0. 

Let's rephrase the above argument in an alternative, coordinate dependent, way. Con- 
sider a unit vector with respect to (11). This vector will have a smaller magnitude with 
respect by (8), as can be immediately seen. So (8) does not increase the magnitudes of 
the unit vectors. Instead, it exponentially decreases their y-component. As a result, it 
will increase distances of any rectifiable curve in the y/transversal direction to x, by an 
exponential factor, which is exactly the statement contained in (20). 

We turn our attention to Lyapunov exponents. We will only be needing their definition 
for the case of uniformly hyperbolic dynamical systems [21], [22], a standard example of 
which is the geodesic fiow on a Riemannian manifold of, generally variable, negative 
sectional curvature. To be slightly more general than that, let ft'.M-^M be a fiow 
on the Riemannian manifold (M, g), whose generating vector field is 

X{t) = j^iMxm^o (22) 

Let Y G TM and let its norm with respect to the Riemannian metric g be indicated 
by ||y|| = {g{Y,Y)}2 . The Lyapunov exponent of the perturbation in the direction of 
Y E T^M along the evolution/trajectory of the fiow ft is defined by 

A,(y)= limilog||4(/tmi (23) 



The Lyapunov exponent measures the asymptotic rate of change of the magnitude of a 
perturbation in the direction of Y . Is is clear that if someone is interested in a stabihty 
analysis of a flow, the most pertinent Lyapunov exponent is the largest positive one in 
some transversal direction. Consider, as a special case, ft to be the geodesic flow on M? 
endowed with (8) that was discussed above. We see from (21) that 

lim ||J(s)|| = e^ (24) 

therefore, the definition (23) gives 

A = 1 (25) 

To summarize: we started from a dynamical system modelled on a Riemannian man- 
ifold M equipped with a metric induced by (11). We assumed that this dynamical 
system's effective statistical description is provided by the Tsallis entropy (1). In turn, 
the Tsallis entropy composition property (6) indicated that it is more suitable to use the 
effective hyperbolic metric (8) on M manifold, instead of the original one (11). It was 
the collective description of the underlying dynamics that dictated this "hyperbolization" 
in producing the effective metric (8) from (11). This "hyperbolization" was concretely im- 
plemented as the "warping" by a convex function, the exponential in the present case, in 
transitioning from (11) to (8). The instabilities of the dynamical system were expressed 
through its positive Lyapunov exponents, initially with respect to (11), the largest of 
them being the most important. When the effective behavior described by the hyperbolic 
metric (8) is taken into account, the largest positive Lyapunov exponent of the underlying 
dynamical system becomes zero, as can be seen in (24) and (25). The reason is that the 
perturbations of the underlying dynamical system and the distances with respect to the 
hyperbolic metric (8) grow at the same rate (25). Hence their relative growth rate is 
linear, as seen in (25), so the Lyapunov exponent of the underlying dynamical system rel- 
ative to (8) is zero. This is the main conclusion of the present work and it is in agreement 
with all known results reached by analyzing particular models [3] (and references therein). 

In case the largest Lyapunov exponent is zero, as in the case discussed above, the in- 
stabilities of the dynamical system grow at a milder than exponential rate, with respect to 
the hyperbolic metric (8), and consequently their asymptotic behavior has to be encoded 
differently, if we want to obtain non-trivial results. One way to quantify the growth of 
such perturbations by modifying the definition of Lyapunov exponents (23) to 

KiY) = hm 1^^^ (26) 

^ ' t^oo hgt ^ ' 



These modified exponents describe perturbations obeying an asymptotic, power-law evo- 
lution t'^'^^'^\ This is essentially the definition adopted in [3], [7] - [9] 

| = A,«' (27) 

for systems described by the Tsallis entropy, in just different notation. Naturally, the 
modified Lyapunov exponents (26) will explicitly depend on the value of the entropic 
parameter q of the dynamical system. Since, as is conjectured, the systems described 
by the Tsallis entropy may possess more than just one value of non-extensive parameter 
q [3], [11], depending on which property of the system is described by them, the modified 
Lyapunov exponents (26) would evidently depend on the value of q determining the 
sensitivity of the system to infinitesimal perturbations, indicated as gsen in [3]- 

The use of a definition like (26) would be of limited interest for applications, if explicit 
constructions/examples did not exist for which it could provide non-trivial information. 
Motivated by the linear (exponential) increase of the geodesic distances in Riemannian 
manifolds of zero (negative resp.) sectional curvature (21), one may be wondering whether 
an intermediate behavior described by (26) is even possible. The answer turns out to be 
affirmative. An example of a quadratic separation of geodesies in a 2-complex endowed 
with a CAT(O) metric was constructeqll in [23]. This quadratic geodesic deviation in a 
CAT(O) space is quite different from the case of CAT(k), k < 0, whose geodesies de- 
viate exponentially from each other, as will be explained in the next Section. Whether 
the property of quadratic or, more generally, polynomial divergence of geodesies encoded 
in (26) is "typical", or even common, for CAT(O) spaces, and whether in such a case 
the definition (26) exhausts all possible geodesic deviation behavior, does not seem to be 
known at this time. More importantly, it is unclear, to us at least, exactly what, if any, 
physical system the CAT(O) space construction of [23] can be used to describe. For these 
reasons, we refer to [23] for the construction itself and its mathematical aspects, aiming 
to re-visit this topic, if physical reasons warrant it in the future. 



3. Geodesic deviation in CAT(k), k < spaces 

We continue with the generalization of the above results to the case of CAT(k), k < 
spaces. The need for using CAT(k), k < spaces, motivated by the composition of the 



^We are grateful to Professor Panos Papasoglou for bringing this work to our attention. 



Tsallis entropy (6), was explained in [12]. The major obstacle in repeating the Rieman- 
nian approach verbatim, is that CAT(k) spaces do not possess a differential structure 
[15], [16], [24], so one has to dispense with statements relying on regularity properties, 
such as ones formulated via vector fields, the geodesic (15) and the Jacobi (16) equations 
etc. The only option left is to employ the triangle inequality, which when combined with 
the CAT(k), k < condition proves to be sufficient for attaining the sought-after goal 
[15], [16], [24]. One should notice, that due to lack of smoothness, most statements in the 
present Section can only be formulated via inequalities, as contrasted to the equalities, 
such as (20), derived in the Riemannian case. As a result, the arguments in this Section, 
applied to the case of CAT(k), k < spaces are, inevitably, synthetic as opposed to the 
analytic ones in the Riemannian case of the previous section. We will, largely, follow [24] 
in the sequel. 

One begins by realizing that a geodesic space {X,d) which is CAT(k), k < is 
necessarily hyperbolic. There are several definitions of hyperbolicity at various levels 
of generality [24] and several equivalences among them [24]. The following definition, 
ascribed to E. Rips, is the most useful for our purposes: a metric space (X, d) is 6- 
hyperbolic, for 6 > 0, if for any three points x,y, z & X, any side of the triangle having 
X, y, z as vertices lies in a (5-neighborhood of the union of the two others. Given this 
definition, the hyperbolicity of the CAT(k), k < space (X, d) follows immediately 
from the definition of the CAT(k) condition. 

Consider a polygon in such (X, d) with n ^^ vertices Xi,a;2, . . .Xn- Let 

+ 1 (28) 



vn 



n 
L2, 



where the square brackets indicate the integer part of their argument. Because (X, d) 
is a geodesic space, every segment has a midpoint. Consider the midpoint Xm and the 
triangle XiXmXn- Since (X, d) is 5- hyperbolic, there is a point y in the union of the 
two other sides XiXm U XmXn such that d{y,z) < 6 where z G x„Xi. Without loss of 
generality, we assume that y G xiXm- Using induction on the m-polygon xi,X2, ■ ■ ■ , Xm, 
we find that the distance of z from the union of all other sides of this m-polygon is 



m—l 



diz, [jx,Xi+A <{p-l)6 (29) 



i=l 



where p G N satisfies 

P>!^4^ (30) 

log 2 

We conclude then, that in (X, d) each side of the rz-polygon is contained in a p6- 

neighborhood of the union of its other sides. 

The second, and last, step is the quantification of the concept of the exponential 
separation of geodesies in {X,d). Consider two segments emanating from xEX toward 
y, 2 G X, respectively. Let two objects move with unit speed, one in each of these 
two arc-length parametrized segments. We are interested in the separation of x and y 
after time t. An obvious way to measure such a separation would be to start at the 
location of y, move back along the segment joining it with a; by a distance t, switch 
segment at x and then continue from x alongside the other segment for an additional 
distance t until reaching z. This is not however, the analogue of a curve joining y 
and z of Section 2. What we want is to measure the "direct" separation between y 
and z, without having to go back close to the intersection x of the two segments. So, 
we want to determine the length of a path from y to z, in the complement of a ball 
Br{x) of radius r > centered at x. Consider a path a;iX2 • • -Xn, n eN, such that 
d{xi,Xi-^-l) < e, i = l,...n — l which lies outside the ball Br{x), with x G x„xi. Let 
JO G N be as in (30). Using the conclusion of the previous paragraph, there exists a point 
w G XiXi+i, i = 1, . . . ,n — 1 such that d{x, w) < pS. Since 

d{x,w) > r — - (31) 



we find 



which gives 



where 



r € 



n>2P-'>c-2^ (33) 



2-(^+i) (34) 



The conclusion that we reach from (33), is that the number of points with a uniform 
distance upper bound e between two consecutive ones, making up the path between y 
and z outside Br{x) increases exponentially with r in the CAT(k), k < space 
(X, rf). This is what we wanted to show. We see that (33) is the analogue of (21) for 
(X, (i). Naturally, the argument leading to (33) is already applicable to the case of a 



Riemannian manifold (M, g) of negative sectional curvature, since in this case (M, g) 
is a CAT(k), k < space. So, the conclusions drawn in the Riemannian case, about the 
effective metric behavior of systems described by the Tsallis entropy, can be extended un- 
altered to the case of any number of interacting systems described by different values of q. 



4. Conclusions and Outlook 

In the present work, we attempted to justify why dynamical systems whose statistical 
behavior is described by the Tsallis entropy, have vanishing largest Lyapunov exponent. 
This was essentially ascribed to employing the effective negative curvature metric (8), 
which is the "hyperbolization" of the Euclidean initially employed metric (11), as was 
pointed out in [12]. Our conclusion is in agreement with all currently known results. 

The process of generalizing this conclusion to the case of CAT(k), k < spaces pre- 
sented in Section 3, is of interest, as it points out to the underlying reason behind such 
behavior. The argument of Section 3 shows that the key in understanding consequences 
of the Tsallis entropy composition property (6) is the concept of hyperbolicity, which was 
also alluded to in [12]. Knowing this, is should not come as a surprise that the argument 
of Section 3 is a small part of a well-known proof of Morse's Lemma, which establishes the 
stability of geodesies in hyperbolic geodesic spaces under quasi-isometries [15], [16], [24]. 
In particular, since Riemannian manifolds of negative sectional curvature (M, g) are a 
subset of CAT(k), k < spaces, the vanishing of the highest Lyapunov exponent of the 
dynamical systems modelled by (M, g) just expresses their underlying hyperbolicity in 
a very compact way. The general framework of this hyperbolicity and its implications for 
systems described by the Tsallis entropy will be examined further in a future work. 
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